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ABSTRACT 

The many similarities between the prompt emission pulses in gamma ray 
bursts (GRBs) and X-ray flares during the fast decay and afterglow phases of 
GRBs suggest a common origin. In the cannonball (CB) model of GRBs, this 
common origin is mass accretion episodes of fall-back matter on a newly born 
compact object. The prompt emission pulses are produced by a bipolar jet of 
highly relativistic plasmoids (CBs) ejected in the early, major episodes of mass 
accretion. As the accretion material is consumed, one may expect the engine's 
activity to weaken. X-ray flares ending the prompt emission and during the 
afterglow phase are produced in such delayed episodes of mass accretion. The 
common engine, environment and radiation mechanisms (inverse Compton scat- 
tering and synchrotron radiation) produce their observed similarities. Flares in 
both long GRBs and short hard gamma ray bursts (SHBs) can also be produced 
by bipolar ejections of CBs following a phase transition in compact objects due 
to loss of angular momentum and/or cooling. Optical flares, however, are mostly 
produced in collisions of CBs with massive stellar winds/ejecta or with density 
bumps along their path. In this paper we show that the master formulae of the 
CB model of GRBs and SHBs, which reproduce very well their prompt emis- 
sion pulses and their smooth afterglows, seem to reproduce also very well the 
lightcurves and spectral evolution of the prominent X-ray and optical flares that 
are well sampled. 

Subject headings: gamma rays: bursts 



1. Introduction 

A flaring activity during the afterglow (AG) phase of a gamma ray burst (GRB) was 
first observed in the late-time AG of GRB 970508 with the Narrow Field Instrument (NFI) 



^ dado@phep3 . technion. ac . il 
Physics Department, Technion, Haifa 32000, Israel 

^ arnon@physics . technion. ac . il 
Physics Department, Technion, Haifa 32000, Israel 



-2- 



aboard the BeppoSAX satellite in the X-ray band (IPiro et al.lll998[). an d with ground based 
telescopes in the optical band (jPian et al.l Il998l : iGalama et al.l |l998aj) . It was interpreted 
in the frame work of the fireb all (FB) model of GRBs as a delayed burst from the central 
GRB engine (jPiro et al.lll998h . Alternatively, in the cannonball (CB) model of GRBs it was 
interpreted as a synchrotron radiation (SR) flare from an encounter of the highly relativis- 
tic jetted ejecta from an underlying superno va (SN) explosion with a density jump in the 
interstellar medium (jPado et al.ll2002l . 120041 ) . Late-time flares were later discovered in the 
broadband AG of several other GRBs that were local ized by the BeppoSAX satellite , most 
notably in GRB 000301C at t ~ 4 days after burst jBerger ei_all I2OO0I: Isagar et al.l I2OO0I ) 



where the flar e was attributed to gravitational microlensing ( Garnavich et al. 200ol ) and 



m 



GRB 030329 (ILipkin et al.ll2004l ) where the flare was interpreted in the framework of the 
FB model as due t o 'refreshed shocks' generated by a late activity of the central engine 
(iGranot et al.ll2003l ). In the CB model, however, these flares were well reproduced by the 
emission of SR from encounters of the jetted ejecta from an underlying SN explosion in a 
star formation region with de nsity jumps within or a t the border of a super bubble created 
by the star formation region (jPado et al.l 120021 . |2004J ). 

Early-time X-ray flares ending the prompt emission phase were also detected by the 
wide field camera (WFC) aboard BeppoSAX in a few GRBs such as GRB 011121. In the FB 
model they were attributed to the onset of the external shock in the circumburst material. 
In the CB model they were interpreted as being due to the last episodes of bipolar CB 
ejections from a shutting off central engine. Shortly after the launch of the Swift satellite 
in November 2004, data collected with its X-ray telescope (XRT) showed that X-ray flares 
are quite common in all the phases of the emission from GRBs. In more than 50% of the 
GRBs observed with the Swift X-ray telescope ( XRT), flares were ob s erved at the end of 



the p rompt emission and/or the early AG phase (IBurrows et al.ll2005l . 120071 : iFalcone et al. 



20071 ). In some cases X-ray flares were observed also at very late times, of the order of 



several days after the prompt emission. Although the information on flares is much more 
sketchy compared to that on the prompt gamma ray pulses, their spectral and temporal 
behaviours show clearly that the X-ray flares during the prompt 7-ray emission follow the 
pattern of the 7-ray pulses, suggesting they are the low energy part of these pulses. The 
X-ray flares ending the prompt emission and those superimposed on the early-time afterglow 
have a fast spectral evolution. Their peak intensities decrease with time, and their spectral 
and temporal behaviours are similar to those of the prompt X/7-ray pulses, except that 
they are progressively softer and last longer. In most cases their 7-ray emission probably is 
below the detection sensitivity of the Swift broad alert telescope (BAT). In some cases their 
f iuence in the XRT band exceeded that of the prompt emission in the BAT 15-150 keV band 
(IChincarini et al.lboOSaE I). 
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Late-time {t^lO^ s) X-ray flares, however, seem to exhibit temporal and spectral be- 
haviours that are different from those of most of the early-time flares. Their power-law 
decline is more moderate, they are "achromatic" with a power-law spectrum almost identi- 
cal to that of the late-time AG, and they show very little spectr al evolution (s ee, e.g., the late 



time broad band flares in GRBs 060614 (IMangano et al. 



20071 ) and 081028 (IMareutti et al. 



20091 ). and the hardness ra tio during la t e tini e flares in Swift GRBs reported in the Swift- 



XRT lightcurve repository, Evans et al.l ( 120091 )). 

Flares in the X-ray lightcurve of Swift GRBs were studied phenomenologically by various 



obseryer groups ([Burrows et al. 



2007: Falcone et al 



2007; 



2005 



2007; 



Chincarini et al. 



Kocevski and Butlerl 120071 : iButler and Kocevski 



2007a 



U l2008aE f). Modifications of previously 



sugge sted models and new theoretical n i odels were proposed and discussed by several au - 
thors JProga et allboosl : king et al.l[2005[ IPai et al.ll2006l iFan et al.l[2006l : IPerna et al.ll2006h , 
but none of these proposed models was shown to actually derive the observed hghtcurves 
and spectral evolution of either early-time or late-time flares from underlying physical as- 
sumptions. 

Flares are a natural consequence of th e cannonball (CB) model of GRBs, which was 
motivated by a GRB-microquasar analogy (jPado et al.ll2002l : iDar and De Rujulall2004l ). In 
the CB model, long-duration GRBs and their AGs are produced by bipolar jets of highly 
relativistic plasmoids o f ordinary matter ejected in acc retion episodes on the newly formed 
compact stellar object (|Shaviv and Dai 



plosions (jPar et al.l Il992l : iDar and Plaga 



995 



parlll998l ) in core-collapse supernova (SN) ex- 



19991 ). It is hypothesized that an accretion disk or 



a torus is produced around the newly formed compact object, either by stellar material orig- 
inally close to the surface of the imploding core and left behind by the explosion- generating 



outgo ing shock, or by more distant st ellar matter falling back after its passage (jPe Rujula 



19871 ). As observed in microquasars (IMirabel and RodriquezI Il999l ). each time part of the 



accretion disk falls abruptly onto the compact object, two jets of cannonballs (CBs) made 
of ordinary-matter plasma are emitted with large bulk-motion Lorentz factors in opposite 
directions along the rotation axis, wherefrom matter has already fallen back onto the com- 
pact object due to lack of rotational support. The entire radiation emitted from a GRB 
is produced by the interaction of the jet of CBs with the environmemt along its path, as 
illustrated in Fig. [TJ The prompt 7-ray and X-ray emission is dominated by inverse Compton 
scattering (ICS) of photons of the SN light filling the cavity produced by the pre-supernova 
wind/ejecta blown from the progenitor star long before the GRB. The CBs' electrons Comp- 
ton up-scatter this 'glory' light into a narrow conical beam of 7 rays along the CBs' direction 
of motion. Each CB produces a single GRB pulse. An X-ray 'flare' coincident in time with 
a prompt 7-ray pulse is simply its low-energy part. The natural explanation of flares ending 
the prompt emission and during the early time afterglow is the same: ICS of glory photons 
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by the electrons of CBs ejected in late accretion episodes of fall-back matter on the newly 
formed central object. Early-time X-ray flares without an accompanying detectable 7-ray 
emission are usually IC flares (ICFs) produced by CBs with relatively smaller Lorentz factors, 
due to weakening activity of the engine: As the accretion material is consumed, one may 
expect the 'engine' to have a few progressively-weakening dying pangs. Like the lightcurves 
of the prompt GRB pulses, the lightcurves of ICFs exhibit a rapid softening during their 
fast decline phase (se e, e.g., the XRT hardnes s ratio reported for Swift GRBs in the Swift 



lightcurve repository, lEvans et al.l ( 120071 . 120091 )). 



In the CB model, each IC flare (ICF) is followed by the emission of SR flare (SRF) from 
the encounter of the CB with the wind/ejecta, which was blown from the progenitor star 
long before the GRB (see Fig. [1]). Because of time-aberration in the observer frame, these 
SRFs appear to have only short lag-times relative to the corresponding ICFs. Below their 
peak energy, the spectral behaviour of the ICS pulses/flares is roughly ~ while the 
SR emission from fast cooling electrons has typically ~ v^^-^ . Thus, while the prompt 
keV-MeV pulses/flares are dominated by ICS of glory light, the 'prompt' optical emission is 
dominated by SR. As the glory extends into the wind, often the SR emission begins before 
the ICs pulse/flare has ended. 

The initial expansion of the CBs and the slow-down of the le ading ones by the circum- 



burst matter may merge most of them into a single leading CB (jPar and De Rujulal 12004 



Dado et al.ll2009al ) during the afterglow phase. Its collimated beam of the prompt gamma 
rays ionizes the matter in front of it. The ions continuously impinging on a CB with a 
relative Lorentz factor 7(t), where 7(t) is the bulk motion Lorentz factor of the CB, generate 
within it an equipartition turbulent magnetic field. The intercepted electrons are isotropized 
and Fermi accelerated by these fields and emit isotropic synchrotron radiation in the CB's 
rest frame, which is Doppler boosted and beamed relativistically into a narrow cone with a 
typical opening angle ~l/7(t) in the observer's rest frame. Late SRFs are produced mainly 
when the CBs encounter winds or density bumps along their path first from the progenitor 
star and later in the interstellar medium (ISM). The lightcurve of these flares depends on 
the unknown density profile of the encountered wind/density bump that cannot be predicted 
a-priori. But, both the early-time and the late-time SRFs have a typical SR spectrum and 
a weak spectral evolution that are quite different from those of the accretion induced ICFs 
and can be used to identify their origin - late ejection episodes from the central engine or 
density bumps. 

In the CB model, short hard bursts (SHBs) are also produced by bipolar jets of plas- 
moids ej ected in mergers of c ompact objects in close binary systems such as neutron stars 
merger (jCoodman et al. 1987 ) or in mass accretion episodes on compact objects in close bi- 
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nary systems, or in phase transit ions in compact stars (neutron stars, hyper stars and strange 
quark stars) ( iDado et aLll2009b| ). Bipolar ejections in late accretion episodes or phase tran- 
sitions after cooling and loss of angular momentum probably produce the observed ICFs in 
SHBs, and SR radiation from encounters with winds/density bumps produces SRFs. 

Flares were routinely included in th e CB model description of the afterglows of GRBs 
and SHBs (IDado et al.ll2002l . l2004l . l2009al Jbl). They were calculated from the master formulae 
of the model, which describe well the prompt ICS emission and the emission of SR at all 
times. It was shown that ICS explains successfully both the prompt keV-MeV pulses and 
the X-ray flares ending the prompt emission, while the SR emitted in the collision of the jet 
of CBs with the wind/ejecta from the progenitor star exp lains well the pronapt optical flares 
meas ured with robotic telescopes in very bright GRBs (IDado et al.l l2009at iDado and Dar 
20081 ) ■ However, attention was focused there on the general properties of the prompt emission 
and the smooth afterglow, rather than on flares during the afterglow phase. Moreover, in 
many ORB afterglows, flares are weak, or are blended, or are not well sampled, and their 
properties could neither be inferred reliably from the AG lightcurve nor used reliably to 
test theoretical models. The situation concerning prominent X-ray flares observed by the 
Swift XRT is different. They are well resolved and their spectral properties are much better 
measured. In this paper we focus our attention on GRB X-ray and optical lightcurves with 
prominent flares. In particular, we compare the CB model lightcurves and their spectral 
evolution with representatives set of X-ray and optical lightcurves measured with the Swift 
XRT, and with ground-based robotic telescopes and Swift UVO in space, respectively, which 
have prominent flares that are well sampled. Such a comparison provides stringent tests of 
both the CB model and its interpretation of the origin of prompt emission pulses in GRBs, 
their afterglows and the early and late time flares in their lightcurves. We show that the 
CB model correctly predicts their main observed properties, and reproduces well their entire 
lightcurves and spectral properties. For completeness, we first summarize the relevant master 
formula of the CB model and their simplified forms that we later use in our analysis of the 
X-ray and optical lightcurve of GRBs. 



2. ICS flares 



In this section we summarize the master fo rmula of the CB mo del for the pulse shape 
and spectral evolution of ICS pulses/flares (see iDado et al.l (l2009al ) and references therein 
for their derivation). Let t denote the time in the observer frame after the beginning of 
a flare. The light-curve of a flare, produced by the electrons in the CBs by ICS of thin 
bremsstrahlung photons filling the cavity formed by a wind blown by the progenitor star 
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long before the GRB, is generally well approximated by (jPado et al.ll2009al ): 



(fN^ tVAt^ E/E^{t) -E/E,{0) p(Ef) (I) 

where A is a constant that depends on the CB's baryon number, Lorentz and Doppler factors, 
and on the density of the glory hght and on the redshift and distance of the GRB, and Ep{t), 
the peak energy of E(PN^/dEdt at time t, is given roughly by: 

E,{t) ^ E,iO) , (2) 

with tp being the time (after the beginning of the flare) when the ICS photon count-rate 
reaches its peak value. For E <^Ep, it satisfies Ep(tp) = Ep where Ep is the peak energy of 
the time-integrated spectrum of the flare. Thus, in the CB model, each ICS pulse in the 
GRB lightcurve is described by four parameters. A, At{E), Ep{0) and the beginning time of 
the pulse when t is taken to be 0. 

Eq. dl]) with Ep given by Eq. ([2]) describes well the shape and the spectral evolution of 
GRB pulses and of early-time X-ray flares. In particular, it correctly describes their rapid 
spectral sof t ening d uring their fast decline as demonstrated in Figs. Wp,f, [3]d and [Tli, and in 



Dado et al.l (j2008bl ) for many more cases. 



If absorption in the CB is dominated by free-free transitions, then At{E) oc-E~°'^, and 
for E<^Ep the lightcurve of an ICF is approximately a function of Et"^ (the law'), with 
a peak at t = At, a full width at half maximum, FWHM ^2 At and a rise time from half 
peak value to peak value, RTf» 0.30 FWHM independent of E. Note that the approximate 
E t^ law makes the fast decline sensitive only to the product Ept^ and not to their individual 
values. This degeneracy in the pulse shape can be removed by inferring Ep from the broad 
band spectrum of the ICF. 

The late-time decay of the energy flux of the prompt emission pulses and ICFs in an 
energy band [El, E2], which follows from Eq. ([T]), is given approximately by: 

E ^{E, t)dE^A Ml^ le-'^'/EAt) _ ^-E2iE,ii)y 

dtdE^ ' ' t2 L J V ; 

Thus, for the Swift XRT lightcurves where El = 0.3 keV and E2 = 10 keV, as long as 
Ep{t) ^£"2 >-E'l, the energy flux in an ICS pulse/flare decays like until it is taken over 
by the SR afterglow. When Ep{t) ^ El but Ep{t)1iE2 the energy flux decays like t""^, and 
when Ep{t)^El the energy flux decays like t"^e~^* /2£^ptp^ 

For a single ICF beginning at t=0, which is superimposed on a smooth SR afterglow, 
the effective photon spectral index V{E,t) = dlog{dN / dE) / dlogE is given approximately by 



flDado et al.ll2008b[ ): 



T{E, t)^[l-(3g + E/Epit)]e{tAG - t) 0(t) + TsR e(t - tAc) , (4) 

where Iag is the time when the SR afterglow takes over the ICS emission, /3g ~ for a glory 
with thin thermal bremsstrahlung spectrum, and Tsr is the best fit photon spectral index 
of the smooth X-ray AG. 

All the above properties are clear fingerprints of ICFs produced by highly relativistic CBs 
fired in mass accretion episodes on the newly formed compact central object. Moreover, the 
observations of X-ray flares indicate that their widths are proportional to their emission time 
tf after the GRB trigger. In the CB model, this implies that their peak energy, equivalent 
isotropic gamma ray energy and peak luminosioty decrease like EpOctJ^ , -EjsoOC [(1 + z)/tf]'^ 
and Lp(x[{l + z)/tf]^ . 



3. The SR lightcurve 



In the CB model, the li ghtcurve in the observer frame of the SR emitted by a CB is 
given by the master formula (jPado et al.ll2002l . l2009al ): 



F,[t] (xA{u,t) 



n{t) R{tf l{t)H{tY 



n-i/2 r 



1 



{l-p)/2 



(5) 



where t is the time after the ejection of the CB, R{t) is its radius, n is the density along 
its trajectory, A{v,t) is the attenuation of radiation along the line of sight to it, //^(t) is the 
typical frequency in the observer frame of the SR emitted by the electrons that are swept 
into it at time t with a relative Lorentz factor 7(t), 

n^/^ -fitf 5{t) 



Mt) oc 



l + z 



(6) 



and p~2.2 is the spectral index of the Fermi accelerated electrons in it. 



3.1. The early time SRFs 



The SR radiation that is emitted from the encounter of a CB with the wind/ejecta of 
the progenitor star, with a density profile n{r) cx: e~"/^''""^'"-'/(r — r,,,)^ for r >r,„ an d n(r) =0 
for r<r^, follows from Eq. and is given approximately by lDado et al.l (l2009al ): 



- 8 - 



where t = T — T^ with T being the time after trigger and the time at the CB-wind 
encounter, tf.r^p is the typical slow-down time of the fast CB expansion, /3 = F — 1, and 
the exponent describes the decreasing attenuation of the emitted radiation when the CB 
penetrates the wind and/or the initial rise in the wind density due to an exponential cutoff 
of the wind ejection. 

Note that fo r ^ tl^^ the asymptotic decline of an SRF is a simple power law 
jPado et aljboosh . 

ocr^z/^r+S (8) 

while that of an ICF is, 

F,[t] cx E-^^ e-^/^-W ~ E-^^ ^-Eit-+tl)/2E,tl 

Thus, ICFs and SRFs may be distinguished by their different tempo-spectral evolution. 

In the X-ray band, early-time ICFs are usually much brighter than their following SRFs. 
But, due to their rapid late-time decay, occasionally the^t"*" tail (Eq. [8]) of the SRF, which 
follows an ICF, can be seen befor e the plateau /shallo w decay phase of the early-time X-ray 
AG takes over (see, e.g.. Fig. 3 in iDado et al.l (l2009al )). 



The early time SR has usually Pox ~ 0.75 and Fy that decreases strongly between 
the optical and the X-ray band. Consequently, although ICS dominates the prompt X-ray 
emission, in the optical band SRFs are usually much brighter than their preceding ICFs, 
which typically have /3~0 around their peak time. Consequently, optical flares are usually 
SRFs. The evolution of the effective spectral index during an optical SRF that follows from 
Eq. ([5]) has the simple form, 

dlogFy[t] 



l/2 + (/3x-l/2) 



(10) 



dlog V ' ■ - ' ' y 

In early-time optical SRFs, Vh that is initially below the optical band can cross above it as n 
of the wind/ejecta increases and then cross back as the density decreases. Such a change in 
/3o, from [3o ^ Px towards /3o ~ 1/2 and back to [jgr^Bx, h as been observed in some early 
time optical flares, e.g., in GRB 071031 kruhler et al.lboOQh . 



3.2. The late-time SR afterglow 



When the merged CBs coast through the constant-density ISM, their SR in the X- 
ray band is well above the cooling frequency o f the Fermi accele rated electrons and their 
unabsorbed AG as given by Eq. ([5]) reduces to (jPado et al. 2009a ): 



1 y-Px 



(11) 
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where 5 = 1/7 {l — Pcos6) is their Doppler factor with 6 being the angle between the hne of 
sight to the CBs and their direction of motion. For 7^ ^ 1 and 6"^ <^ 1, 6 ^ 2'-/ / (1+7^ 0"^) to 
an excellent approximation. In the CB model the canonical value of the spectral index well 
above the bend frequency Ub has the value /3x ~ 1.1 . For a CB of a baryon number N^, a 
radius R and an initial Lorentz factor 70, relativistic energy-moment um conservation yi elds 
the deceleration law of the CB in an ISM with a constant density n ( iDado et al.ll2009al ): 



7(t) 



7o 



(12) 



with to = (1 + 2;) iVg/S cn TT i?^ 7o . As can be seen from Eq. f|T2l) . 7 and 6 change little as 
l ong as t <C th = [1 + 796*^]^ to, and Eq. (fTTl) yields the 'plateau' phase of 'canonical AGs' 
(INousek et al.ll2006l ). For t^th, 7 and S decrease like . The transition 70 -^7o (^Ao)"^^"^ 



around tf, induces a bend, the so called 'jet break', in the synchrotron AG from a plateau to 
an asymptotic power-law decay. 



FisAi[i^,t] oc r 



-1/2 ^-fc 



(13) 



Thus, the shape of the entire lightcurve of the SR afterglow from a CB that enters the 
constant density ISM depends only on three parameters, the product 70 9, the deceleration 
parameter (or the break time th) and the spectral index Px- The post break decline is 
given by the simple power- law (Eq. [13]) independent of the values of 70 9 and t^. In cases 
where tb is earlier than the beginning of the XRT observations or is hidden under the prompt 
emission, the entire observed lightcurve of the AG has this asymptotic power-law behaviour 
JPado et al.ll2008ah . 



3.3. Late-time SRFs 

The lightcurve of late-time SRFs strongly depend on the density profile of the density 
bumps. For a wind-like density jump ?T,oce~"™/*^''"^™-'/(r — r^)^ beyond r = r^, the lightcurve 
is given by Eq. ([7]). Such a density profile is expected for the boundaries of star formation 
regions where GRBs usually take place. At late times both the optical and the X-ray bands 
are above the bend frequency, /So ~ and the asymptotic decline of the corresponding 
late-time SRF is given by: 

F,[t](xt-^u-^-\ (14) 

where t = T — is the observer-time after the CB has reached r^. At late times, the bend 
frequency is well below the optical band and remains so during the crossing of density bumps. 
Thus, no detectable spectral variation is expected in late-time SRFs. 
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3.4. Correlations 

Because of the large bulk motion Lorentz factors of the jets of CBs, Doppler boosting, 
relativistic beaming and time aberration yield strong dependence of their observed radia- 
tions on 7 and 6. This dependence dominates the flare observables and can be used to 
correlate triplets of independent flare observables without knowing the exact values of 7 and 

5. Moreover, many independent observables depend on the same combinations of 7 and 

6, which results in pair correlations. Finally, due to selection effects, various observables 
depend strongly only on the Lorentz factor or on the Doppler factor, which also yields pair- 
correlations. These correlations between various radiations, between flare observa bles and 



between flare and GRB observables are discussed in detail in iDado and Pari (120101 ). 



4. Case Studies 

Although in the CB model early-time optical SRFs follow the X-ray ICFs, they have a 
different origin and a different tempo-spectral evolution. Thus, we shall compare separately 
the observational data on X-ray flares and optical flares and the CB model predictions. 



4.1. X-ray flares 



In the CB model, the X-ray emission in GRBs, XRFs and SHBs during the prompt 
emission phase is a sum of X-ray flares that are part of their pror npt gamma r ay pul ses. 
That has been shown repeatedly in CB model publications (see, e.g., iDado et al.l (l2009al Jbl)) 
and independently by many other authors. We have also shown that the fast decline phase 
of their prompt emission with its rapid spectral softening is just the tail of the last prompt 
emission pulses, because the exponential factor in Eq. ([1]) suppresses very fast the relative 
contribution of the earlier pulses by the time the data sample the later pulses or flares. This 
i s demonstrated in F igs. 2a,b,c where we compare between the 'canonical X-ray lightcurve' 
( Nousek et al.l 120061) of GRB 060729 and the evolution of its spectral index as inferred by 
Zhang et al.l (120071 ) and the C B model de s criptions of th ese lightcurves. This has already 
been demonstrated in detail in iDado et al.l (j2008bl l2009al ) for many oth er GRBs and XRFs . 
It is also the case for SHBs as shown in Figs. 2e,2f for SHB 050724 and in lPado et aD J2009bl ) 
for many other SHBs. GRBs with a non-canonical X-ray AG, are simply GRBs where the 
emission of SR begins to dominate the X-ray lightcurve before the fast decline phase of the 
prompt emission. Such a case is shown in Fig. 3d where we compare the XR T lightcurve 
of GRB 060418 and its CB model description. Many other cases are shown in iDado et al. 
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fl2008bl . l2009al Jbh. 



In order to have a stringent test of the CB model interpretation of X-ray flares during 
the afterglow phase of GRBs, we have selected a sample of 14 GR Bs with X-ray lightcurves 
that are reported in the Swift/XRT GRB lightcurve repository (lEvans et al.l 120071 . 120091 ) 
and have well sampled prominent flares during their AG phase. This sample includes GRBs 
050502B, 050916, 060526, 060929, 070704, 080506, 080607, 080810, 080906, 081102, 090417B, 
090621A, 090709A, and the SHB 050724. For all these GRBs we have fitted the entire XRT 
lightcurve with the master formulae of the CB model. In order to minimize the number of 
adjustable parameters ir i the theoretical li ghtcurves, we adopted the standard simplifying 
CB model assumptions ( jPado et al.ll2009a( ): The burst environment is a cavity full of thin 
bremsstrahlung optical photons (glory light) enclosed within a wind/ejecta that have a den- 
sity profile noc e""""/^^"^'""-*/(r— r^;)^ beyond until the density is taken over by the constant 
density of the ISM. The CBs were taken to be well separated during the prompt emission 
phase and to be well represented by a single effective CB during the afterglow phase. The 
latest one or two observed pulses/flares in the prompt emission were assumed to dominate 
its fast decay. This fast ICF decay is overtaken by the tail of the SR emission from the 
encounter of the CBs with the wind/ejecta as given by Eq. ([7]), or by the plateau phase of 
the SR afterglow emitted from the decelerating CB in the constant-density ISM. The smooth 
afterglow was calculated from Eqs. f lTTj) . ( !T2|) . using best fit values of the normalization, 70 9, 
to, and p/2 =(3x with (3x wit hin the error range reported for T = f3x+l in the Swift X-ray 
repository (lEvans et al.ll2009l ). 



Prominent X-ray flares during the afterglow with a rapid spectral softening during their 
decline phase were assumed to be ICS pulses/flares. Such flares were superimposed on the 
CB model smooth SR afterglow. Flares with a constant hardness ratio similar to that of 
the smooth SR afterglow were assumed to be SRFs. They were generated by introducing 



density bumps with a windy profile into the master formula of the SR lightcurve (jPado et al. 
2009ah . 



The fitted parameters of the CB model descriptions of the 14 X-ray light curves of the 
above GRBs are listed in Tables \Tf2[ When only the tail of the SRF was visible we used 
its parameter- free asymptotic form, Eq. ([H]). Because of the use of simplifying assumptions, 
the values of the parameters may be effective values and not true values. Therefore, we 
refer to the CB model fits as 'descriptions' rather than as best fit predictions. In order to 
avoid repetitions and an excessively long section, we limit our detailed discussion to three 
representative cases, GRB 050502B representing long soft GRBs, SHB 050724 representing 
short hard bursts, and the recent peculiar GRB 090709A with a suspected 8 s periodicity in 
its prompt emission. 



-VI- 



GRB 050502B was studied in detail by lFalcone et al.l (|2006[ ). The XRT began taking data 
63 s after the BAT trigger and followed its X-ray lightcurve until 10.6 days after burst. The 
measured lightcurve in the 0.3-10 keV band is shown in Figure 2a. Following an initial low- 
flux level, the XRT detected a giant X-ray flare, which began at 345 ± 30 s, reached a peak 
value around 770 s, with intensity more than 500 times larger than that of the underlying 
afterglow. The fiuence of the flare, (1.0±0.05) x 10^ erg cm~^ in the 0.2-10.0 keV energy band, 
exceeded that of the prompt emission measured with the Swift BAT in the 15-350 keV energy 
band. After several hours, two weaker flares in the X-ray emission occurred consecutively 
beyond which the decay of the X-ray lightcurve became steeper. Except for the giant flare, 
the spectrum of the afterglow and the two late flares was well fit by a power-law with a 
photon spectral index 1=1.945 (+0.077, -0.100) jEvans et al.ll2007l . boooh . The spectrum of 
the flare was well fit with a power-law with an exponential cutoff; however, due to its non 
detection by the BAT, its value could not be well determined from the spectrum measured 
only by the XRT. The photon spectral index of the giant flare before peak-time was F^l, 
much harder than that measured before and after the flare. During the fast decline phase of 
the flare, a rapid spectral softening took place and the spectral index increased rapidly to a 
value well above the constant value of the smooth afterglow. 

CB model fits to the X-ray lightcurve of GRB 050502B and the evolution of its spectral index 
are shown in Figs. 3a,b. They show that the observed lightcurve and spectral evolution of 
the giant X-ray flare are well described by the master formula (Eq. ([1])) of the CB model 
for an ICF. As can be seen in Figs. |3^,b, the fast decay and spectral softening stopped 
simultaneously when the AG was taken over by the smooth SR from the decelerating CB 
in an ISM of a constant density. The late-time afterglow shows two flares superimposed on 
the smooth AG, similar to those observed in many long GRBs. Their spectrum, which is 
similar to the smooth AG, suggests that they were produced by enhancement of the emitted 
SR when the CB encountered density bumps in its voyage through the host's ISM. Such late 
flares with a typical SR spectrum and little spectral evolution produced by density bumps 
along the CB trajectory in the ISM add no specific information on the origin of GRBs. 



SHB 050724 at redshift z = . 257 w as studied in detail by ICampana et al.l ( 120061 ). citet- 



Grupe2006 and iMalesani et all ( 120071 ). The Swift BAT triggered on the burst at 12:34:09 
UT on July 24, 2005. The burst had T9o=3.0±1.0 s, but most of the energy of the initial 
SHB was released in a hard spike with a duration of 0.25 s. The bulk of the burst energy 
was not emitted in the short initial spike but in an extended soft emission component that 
lasted ~ 150 s. Swift XRT began observing the afterglow 74 s after the BAT trigger. The 
Chandra X-ray observatory performed two observations, two days and about three weeks 
after the burst. The complete X-ray lightcurve is shown in Fig. 2e. It has a rapid decay with 
a fast spectral softening ending with a sharp transition to a shallower decay with a much 
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harder power-law spectrum, r = 1.79±0.12 (Swift repository, lEvans et al.l (l2009l )). The AG 
steepens gradually into a late power-law decay. A large flare superimposed on the canonical 
lightcurve occurred around 50 ks after burst with a fiuence of ~7% of that of the prompt 



burst. The flare has been detected also i n the optical and NIR bands, e.g., iMalesani et al. 



( I2OO7I ). Spectral analysis of the XRT data (ICampana et al.ll2006l ) showed no evolution during 
the afterglow phase, including the large late flare. Spectral ar ialysis of the Chan dra observa- 
tions from the fading tail of this flare confirmed this result (IGrupe et al.ll2006|). The burst 



took place 2.5 kpc (in projection) from the center of an elliptical host galaxy ( IMalesani et al. 
2007h . 



CB model fits to the X-ray lightcurve of SHB 050724 and the evolution of its spectral index 
are shown in Figs. 2e,f. The initial fast decay with a rapid spectral softening is described 
in the CB model by the tail of an ICF. As can be seen in Figs. 2b, 2c, the fast decay and 
spectral softening stopped simultaneously when the AG was taken over by the SR from the 
decelerating CB in an ISM of a constant density. As expected, the late-time SR afterglow 
is similar in shape to the SR afterglow of long GRBs. Also the late-time flare superimposed 
on the canonical lightcurve is similar to those observed in many long GRBs. Its spectrum, 
which is similar to the smooth AG, suggests that it was produced by enhancement of the 
emitted SR when the CB encountered a density bump in its voyage through the host's 
ISM. The Chandra data show that the canonical AG continued to decay after the flare with 
the same slope and the same spectral index, = 0.79 ±0.15, as that of the AG before 
the SRFs. As shown in Fig. 2e, the complete XRT lightcurve is well described by the CB 
model. Moreover, the CB model relation for the asymptotic decline of the smooth AG is 
well satisfied: The temporal behaviour of the smooth AG was best fit with p=1.56, implying 
an unabsorbed spectral index, /3j>c=p/2=0.78, in agreement with that inferred from the XRT 
and Chandra observations. The observed elliptical host galaxy of SHB 050724 was argued 
to provide strong support for a neutron star merger origin of this SHB. But, it was pointed 
out that neutron star mergers do not produce the late accretion episodes neede d to power a 



late c entral activity that could produce the large flare around 50 ks after burst ( IGrupe et al. 



20061 ). In the CB model, a late flare with a typical SR spectrum and little spectral evolution 
is produced by density bumps along the CB trajectory in the ISM. Such flares neither rule 
out nor support any specific origin of the SHB. 



GRB 090709A triggered the Swift BAT on July 9, 2009 at 07:38:34 UT (IMorris et al. 



20091 ). Its prompt emission lightcurve within 100 s econds after trigger ap peared to show a 
quasi periodic variation with a period of~8 seconds (IMarkwardt et al.ll2009l). which may have 



also seen in indepen dent measuremen ts of its hghtc urve with Konus-W ind (iGolenetskii et al. 



2009h . INTEGRAL JCotz et al.ll2"o09[ ) . and Suzaku Johno et al.ll2009h . The Swift XRT began 
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its follow-u p observations 68 seco n ds af ter the BAT trigger (IMorris et al.l 120091 ). However, 
analysis by iMirabal and Gotthelj ( 120091 ) of the XRT data during 79-469 seconds after the 
BAT trigger did not reveal any significant periodicity. From deep optical observations that 
did not detect its host galaxy, it was concluded th at GRB 090709A took pl ace either in 
the Milky Way or at a redshift between 8 and 10 (jCastro-Tirado et al.ll2009l ). The large 
redshift and the fluence of 9.1 x 10~^ergcm~^ measured with Konus-Wind in the 20 keV-3 
MeV energy range imply an isotropic equivalent gamma-ray energy between 8.6 x 10^^ and 
1.1 X 10^^ ergcm"^, which makes GRB 090709A more luminous than any GRB with known 
redshift. 

Figs. [6^,f present a comparison between the XRT lightcurve of GRB 090709A and its 
CB model description assuming it was an ordinary GRB. The early time flare around t=90 
seconds probably is the X-ray part of the prompt emission pulse with a peak around 87 
seconds in the Swift BAT lightcurve that is slightly delayed relative to the gamma-ray peak, 
as expected from the Et"^ law of the CB model. Its fast decline and rapid softening are 
those expected from an IGF. They are taken over by SR around 150 seconds after the BAT 
trigger. The next 3 peaks probably are SR peaks as suggested by their shape and their 
hardness ratio that is roughly the same as that of the SR afterglow. The best fit CB model 
lightcurve of the smooth AG yields p=2.16, which, in the CB model, implies x = 1.082, i n 



good agreement with the best fit spectral index reported in the Swift repository (lEvans et al. 
2009h . /3x = r-l =1.081 (+0.076, -0.074). All together, GRB 090709A looks like a normal 
GRB with a normal early time flaring activity that took place at a relatively very large 
redshift, and its X-ray hghtcurve is well reproduced by the CB model. 



4.2. Optical flares 



The CB model predictions for early-time and late-time flares in the optical lightcurves 
of GRBs are compared with observations for a representative set of GRBs in Figs. [TlfTTl For 
all these GRBs we show both their measured optical and X-ray lightcurves (if available) and 
their CB model descriptions with the parameters listed in Tables [TE 

Only in a very few bright GRBs w as the prompt optica l emission resolved into separate 
flares . Two such cases, GRB 080319B JWozniak et al.ll2009h and GRB 071003 JPerlev et al. 
20081 ) are shown in Fig. [71 In most GRBs and XRFs, the prompt optical emission that 
appears as a single extended flare, probably, is a sum of unresolved flares. Cases where 
the prompt optical flare was partially resolved into a s um of flares or wher e there is clear 
evidence for overlap ping flares are, e.g., GRB 080330 (IGuidorzi et al.l 120091 ) . XRF 071031 
kruhler et al.lbood ). GRB 061007 favkoff et al.ll2009h . which are shown in Figs. Eli In 
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most GRBs where the prompt optical emission was detected with robotic telescopes from the 
ground or from space (with the Swift/UVO), the prompt emission appears like a single flare, 
probably as a result of either strongly overlapping flares, or insufficient temporal resolution 
due to low statistics, or because the prompt emission did con sist of a single flare . Exam- 



ples of GRBs with a 'single' p rompt optical flare ar e 990123 (lAkerlof et al. 



( Rvkoff et al 



( iGendre et al. 



2004*), 050820A JVestrand et all bood). 081 203A ( 



Kuin et al 



19991 



20091) 



030418 
090102 



20090, 090618 Jli et al.lT2009[ l and 091029 JlaCluvze et al.ll2009[ l. which are 



shown in Figs. IWTTl Examples of brigh t GRBs with well r esolved late-time op tical flares 
include GRB 030329 Jlipkin et al.ll2004l ) and GRB 060206 JWozniak et al.ll2006l ) shown in 
Fig. [nf. However, in order not to inflate this paper, only GRB 080319B, XRF 071031, GRB 
061007 and GRB 990123 are discussed here in detail. 

GRB 080319B at redshift z = 0.937, the brightest GRB observed so far,was simultane- 
ously detected by the Swif t-Burst Alert Telescope (BAT) and the Kon us gamma-ray detector 
aboard the Wind satellite (IRacusin et al.ll2008l : iGolenetskii et al.ll2008l ). The location of GRB 
080319B was fortuitously only 10° away from GRB 080319A, which was detected by Swift 
less than 30 minutes earlier, and allowed several wide fleld telescopes to detect the optical 
emission of GRB 080319B instantly. It started after the beginning of t he prompt keV-MeV 
emission and it peak ed 26 s after the Swift trigger at magnitude 1^ = 5.3 ( IRacusin et al.ll2008 



Wozniak et al.ll2009[ ) visible to the naked eye. The extreme brightness of the burst and its 
gamma-ray. X-ray and UVOIR afterglows led to a flurry of follow-up ob servations wi t h a va - 
riety of space- and gr o und-based telescopes, which were summarized by iBloom et al.l ( 120081 ) ; 
Racusin et al.l ( 120081 ) ; IWozniak et al.l ( 120091 ) . Its isotropic equivalent gamma-ray energy re- 
lease was Eiso^l-'i X 10^^ erg, similar to that of GRB 990123. The fast spectral variation of 
its hard X-ray and gamma ray emission was well parametrized with an exponentially cut-off 
power-law with a cut-off energy that was strongly correlated with the peak structure of the 
lightcurve and a low-energy photon spectral index, Fr^ 1, which changed abruptly into r^2.1 
after the fast decay phase of the prompt emission. The optical and gamma-ray ligh tcurves 



during the explosion were not correlated (see, e.g.. Fig. 1 in IRacusin et al.l (120081 )): The 
onset of the optical emission lagged behind the gamma ray emission by several seconds and 
decayed more slowly at the end of the prompt emission. The typical time scales of their 
temporal variability were entirely different. The extremely bright optical emission could not 
be reconciled with a single emission mechanism - extrapolating the gamma-ray spectrum to 
the optical band underestimates the optical flux by more than 4 orders of magnitude. Their 
spectra were also quite different. Contrary to flreball (FB) model expectations, the X-ray and 
UVO afterglow lightcurves were also chromatic, with no obvious 'jet breaks' and with spec- 
tral and temporal power-l aw decays th a t did n ot satisfy the c l osure re lations expecte d in th e 



FB model (see, however, iBloom et all (120081 ): IRacusin et al.l (l2008l ): IWozniak et al.l (120091 ): 
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Kumar et al.l ( 120081 ) for attempts to reconcile the observations with the fireball model). 



The prompt 7-ray and hard X-ray e r nissio n in GRB 0803 19B is composed of many 
narrow peaks (see Fig. 1 in iRacusin et al.l (120081 ) ). most of which are not well r esolved. Its 
0.3-10 keV X-ray lightcurve measured with the Swift XRT ( IRacusin et al.l 120081 ) and its CB 
model description assuming a constant ISM density until around t~4 x 10^ s, presumably 
when the CB escaped the star formation region into the halo of the host galaxy, are shown 
in Fig. 7a. The best fit parameters are listed in Table HJ The late time temporal decay 
of the X-ray AG is well described by a power-law with ax = 1.54±0.04, except around 40 
ks, where the lightcurve is poorly sampled. This value of ax satisfies well the CB model 
closure relati on ay=/jy-t-l/2=1 .53! =b0.064. As expected for GRBs with large measured Ep, 
Eiso and Lp (jPado et al.ll2008bl ). no AG break was observed in the XRT lightcurve before 



6 X 10^ s. The wiggling of the measured lightcurve around a power-law decay is probably 
due to variations in the ISM density along the CB trajectory, which we have not tried to 
parametrize. A late-time SRF with a typical late-time decays ii~tfY^'^ probably was 
observed around 6 x 10^ s. 

In Fig. 7b we compare th e measured -R-band (and V band renormalized to the R band) 
lightcurve of GRB 080319B (IRacusin et al.l 120081 ) and its CB model description assuming 
that the initially expanding 3 CBs merged into a single CB by the end of the prompt ICS 
emission of gamma-rays and hard X-rays around 300 s (observer time), which decelerates in 
roughly a constant density ISM. The afterglow parameters are listed in Table |3l The 'missing 
jet break' probably is hidden under the prompt emiss ion. Shown also is th e contribution to 
the i?-band afterglow from an SN akin to SN1998bw ( IGalama et al.lll998bl ) displaced to the 
GRB site. 

The early-time optical emission that was resolved into three prominent peaks and its CB 
model description in terms of 3 SR peaks, each one described by Eq. ([7|) with the parameters 
listed in Table |3l are compared in Fig. 7c. The decay of the prompt emission can also be 
reproduced assuming a single CB crossing 3 wind shells, which were ejected by the progenitor 
star long before its SN explosion, rather than 3 CBs crossing a continuous pre-supernova 
wind blown by the SN progenitor. 

In Fig. 7d, the observed evoluti on of the spectral ind ex in the 15-150 keV band during the 
prompt emission and the afterglow (IRacusin et al.ll2008l ) are compared with that predicted by 
the CB model. The predicted sharp transition from the prompt emission, which is dominated 
by ICS of thin bremmstrahlung (Ffal), to SR (ra.Ri2.1), which dominates the emission once 
the CBs encounter the progenitor's wind/ejecta, is clearly observed. 



XRF 071031 was studied in detail by Kruhler et al. (2009). The Swift/XRT began follow- 
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up observations of X RF 071031 1 03 s after the burst. The early XRT hghtcurve (XRT 
hghtcurve repository, Evans et aLl (j2009[ ) is dominated by bright flares at around 120 s, 150 
s, 200 s, 250 s and 450 s. The late X-ray data exhibit re-brightenings at 5.5 ks, 20 ks and 
55 ks superimposed on a smooth AG. The complete XRT hghtcurve was reproduced with 
the CB model, assuming 7 early-time ICFs plus two late time SRFs superimposed on a 
smooth SR afterglow. Their parameters are listed in Table [2l As shown in Figs. [8t,d the CB 
model describes well the very complex XRT hghtcurve (x^/rfo/ = 431/426). The values of 
the photon spectral index of the late time afterglow, rx=1.86 (-f0.14, -0.15), and the index 
of the late time power-law decay, ax = 1.4 ± 0.1, which w'ere inf erred from the Swift/XRT 
observations (Swift XRT hghtcurve repository. lEvans et al.l ( 120091 ) . satisfy well the CB model 
closure relation, ax = rx-l/2. 

The emission in the optical ban d was detected and followed up from the ground in 
automated observations by GROND (IKruhler et al.l l2009l ) with good temporal resolution, 
which began 225 s after burst and lasted nearly 7 h. The white hghtcurve obtained by 
adding the various optical bands, in order to increase sensitivity, shows a broad peak with 
clear deviations from a power-law rise and decay. Fig. [8^ shows the CB model description of 
this hghtcurve in terms of strongly overlapping 4 SRFs following their preceding ICFs in the 
X-ray hghtcurve. The parameters used in the CB model description are list ed in Table El 



Comp arison between the optical spectral index hghtcurve that was inferred by lKruhler et al. 
(120091 ) and the CB model prediction as given by Eq. (|T0|) is shown in Fig. |8f. 



GRB 061007 at a redshift z=1.26 (jJakobsson et al.ll2006l ). which was detected by Swift, 
had a T90 duration of 75.3 s in the Swift BAT 15-350 keV band, and consisted of three 



large peaks with a long faint tail (IMarkwardt et al.ll2006t ISakamoto et al.l 120081 ). The Swift 
XRT began observations 80 s after the start of the burst. The optical observations began 
27 s after the start of the burst by the ROTSE-IIIa robotic telescope (IRykoff et al.l 120091 ) 
The optical hghtcurve shows a sharp rise, brightening by over a factor of 50 in less than 5 
s, followed by two peaks and a steady power-law decline. The 7-ray and X-ra y lightcurves 



have multiple peaks that are not contemporaneous with the optical peaks (jRykoff et al. 



2OO9I ). followed by a steady decline in the X-rays that tracks the optical decline. The three 
overlapping optical peaks could be the three SRFs following the three large BAT peaks. 
In Fig. int we compa re the 0.3-10 keV S wift XRT hghtcurve of GRB 061007 (Swift XRT 
hghtcurve repository, lEvans et al.l (120091 ) ). and its CB model description as the tail of the 
prompt emission SRFs overtaken around i:=300 s by a standard AG with an early-time break 
hidden under the prompt emission. The best fit parameters are listed in Table [1] Note that 
both the early-time and late-time decays satisfy well the CB model closure relations: a = 
r = 1.884(+0.023, -0.023) for the early-time decay and a = r-l/2 = 1.518(+0.087, -0.083) 
for the late-time decay, respectively. In Fig. |9ll we compare the optical hghtcurve of GRB 
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061007 as measured by ROTSE-IIIa flRykoff et al.l l2009l ) and its CB model description in 
terms of three overlapping SRFs and a standard AG taking over around t=400 s. Note 
that the late-time optical AG also satisfies well the CB model prediction, ao = /^o + 1/2 = 
/3x + l/2 = 1.518(+0.087, -0.083). 



GRB 990123 at redshift z = 1.600 (IKulkarni et al.lll999l : lAndersen et al.lll999l ) was for a 
long time the brightest known GRB. It was also the first GRB in which an optical emission 
was detected during the prompt 7/X-ray emission. GRB 990123 was detected and localized 
by the Burst And Transient Source Experiment (BATSE) on board the Compton Gamma 
Ray Observatory (CGRO) in the keV-MeV range an d at higher energies by the COMPTEL, 
OSSE and EGRET instruments (jBriggs et al.lll999l ). It was also detected and localized by 



the Gamma Ray Burst Monitor (GBM) aboard the BeppoSAX satellite (IMaiorano et al. 
2OO5I ) . The lightcurves of the prompt emission in the keV-MeV range showed a complex 



structure of at least 9 pulses/fiares (see Fig. The prompt optical emission that was 

detected by the Robotic Optical Transient Search Experiment (ROTSE) at Los Alamos 22 
s after the onset of the burst, brightened and pe aked at magnitude V ^^9, about 50 s after 
the GRB onset, and decayed rapidly with time (lAkerlof et al.lll999l ). The prompt optical 
fiare was not resolved int o separate fiares. It was followed in the UVONIR bands with large 



Holland et al 



ground based telescopes (ICastro-Tirado et al.lll999l : iGalama et al.lll999l : iKulkarni et al.lll999 



2000h Hu bble Space Telescope until it faded 

27.7±0.15, two months after burst (IFruchter et al.ll2000l ). The broad 



Fruchter et al.lll999 
to a magnitude V = 

band 7-ray, X-ray, UVO and NIR lightcurves of GRB 990123 were reanalyzed recently 
within the synchrotron fireball (FB) model by Corsi et al. (2005). Essentially they found 
that the spectral and temporal properties of the prompt optical emission are uncorrelated 
to the 7 and X-ray emission, implying different physical origins, that the optical and X-ray 
afterglow lightcurves are chromatic contrary to expectations, and that their spectral and 
temporal power-law decays do not satisfy the closure relations of the FB model. 

In Fig. [9k we compare the BATSE multi-peak lightcurve of GRB 990123 in the 20-50 
keV channel (jBriggs et al.lll999l ) and its CB model description. The count-rate in the 20-50 
keV energy band was calculated from the integral J F^dE/E using Eq. ([1]) with the best 
fit parameters, which are listed in Table H] for the 9 peaks sugges ted by the multichannel 
BATSE data and by the BeppoSAX data (IMaiorano et al.l l2005l ). As shown in Fig. [Ht^, 
the shape of the peaks and the entire lightcurve are well reproduced by Eq. ([1]). The 2-10 
keV lightcurve of the X-ray afterglow of GRB 990123 that was measured with BeppoSAX 
(IMaiorano et al.l |2005| ) for t < 2.5 days (not shown here) was best fit by the CB model 
with p = 1.79, implyin g I3x = 0.90, consistent with /3x = 0.94 ±0.12 that was inferred by 



Maiorano et al.l (120051 ) from their data. The ob served temporal powe r-law decay index of 
the late-time X-ray afterglow, ax = 1.46 ±0.04 (IMaiorano et al.l 120051 ) . also obeys the CB 
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model relation , ax =/3x + l/2 = 1.44±0.13 . 



In Fig. |9b we compar e the observation s of the optical light curve of GRB 990123 frorn 



onset (Akerlof et al. 19991) until late tim e ( Castro-Tirado et al. 1999; Galama et al 



Kulkarni et al. 



1999 



Fruchter et al. 1999; Holland et al. 



1999 



20001). normalized to the l^-band, 



and its CB model description as given by Eq. (fTT]) with the afterglow parameters 7 6' = 0.24, 
to = 2250 s and p = 1.79 . Due to a gap in the data between 500 s and 15,000 s, the expected 
transition from a circumstellar density profile cx 1/r^ to a constant ISM density was not well 
determined. However, the gradual bending ('jet break') of the optical AG to an asymptotic 
power-law decay, (xt~^'^~^^'^u~^° , is well reproduced with the expected late-time spectral 
index /3o~/3x~l-l- 



5. Summary and Conclusions 

In the CB model, GRBs, XRFs and SHBs and their afterglows are produced by the 
interaction of bipolar jets of highly relativistic plasmoids (cannonballs) of ordinary matter, 
which are ejected in mass accretion episodes on a newly formed compact stellar object, with 
the radiation and matter that they encounter along their path. As observed in microquasars, 
each time part of the accretion disk falls abruptly onto the compact object, two jets of 
cannonballs (CBs) made of ordinary-matter plasma with large bulk-motion Lorentz factors 
are emitted in opposite directions along the rotation axis, wherefrom matter has already 
fallen back onto the compact object due to lack of rotational support. The prompt 7-ray 
and X-ray emission is dominated by inverse Compton scattering (ICS) of photons of the glory 
- a quasi isotropic optical light emitted by the supernova and scattered by the wind/ejecta 
blown from the progenitor star long before the GRB. The CBs' electrons Compton up-scatter 
the glory photons into a narrow conical beam of 7 rays along the CBs' direction of motion. An 
X-ray 'flare' coincident in time with a prompt 7-ray pulse is simply its low-energy part. The 
early-time X-ray flares without a detectable accompanying 7-ray emission are usually ICFs 
produced by CBs with relatively smaller Lorentz factors due to a weakening activity of the 
central engine: As the accretion material is consumed, the 'engine' has a few progressively- 
weakening dying pangs. The hghtcurves of ICFs, like those of the prompt emission pulses, 
exhibit a rapid softening during their fast decline phase. Roughly, the lightcurves are a 
function of the product Et"^ and not of the individual values of the photon energy E and 
the time t after the beginning of the flare. The peak energy, isotropic equivalent energy and 
peak luminosity of the ICFs are correlated like those of the prompt GRB pulses. 

In the CB model, each ICF is followed by SRF from the encounter of the CB with 
the wind/ejecta that was blown from the progenitor star long before the GRB. Because of 
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time-aberration, in the observer frame, these SRFs lag after their preceding ICFs by a short 
time of the order of the ICS pulse duration. Optical flares are usually much wider than their 
corresponding gamma/X-ray pulses and overlap, which makes it difficult to associate the 
early time optical flares with their preceding gamma/X-ray pulses/flares and measure their 
lag-time. Only in single-pulse GRBs that are bright enough to be detected with robotic tele- 
scopes and/or Swift UVO and in very bright GRBs, such as 080319B and 071031, where the 
optical flares were partially resolved with robotic telescopes, could the predicted association 
between early time optical flares following gamma/X-ray pulses/flares be tested. 

Often the fast decay of an X-ray ICF is taken over by SR of X-rays from the CB encounter 
with the wind enclosing the glory light be fore it disappears un der the plateau/shallow decay 



phase of the AG (see, e.g.. Figs. 3a,b,c in iDado et al 



fl2009al )l 



Late-time flares are usually SRFs produced by CB encounters with the bumpy boundary 
of the star formation region, or with density bumps within this region where the GRB took 
place. Their exact proflle is not known apriori but a wind-like proflle seems to be a good 
working hypothesis. 



Unlike the empirical parametriz ations (such as Band fu nction ( iBand et al.lll993[ ). cut 



off power-law, Beuermann function (IBeuermann et al.lll999l ). broken power- law, segmented 



power-law, etc.) used in most of the published standard analyses of GRB lightcurves, which 
have never been properly derived from underlying physical assumptions, the master formulae 
of the CB model were derived in fair approximations from its underlying physical assump- 
tions. As shown in this paper, the lightcurves and spectral evolution of X-ray flares and 
optical flares in GRBs, XRFs and SHBs are well described by the master formulae of the 
cannonball model of GRBs. So far, no new assumptions or modiflcations of these formulae 
were needed when applied to well sampled flares in the GRB lightcurves. Probably, in the 
future, when much more reflned data and better sampled lightcurves of GRBs and their AGs 
will become available, the CB model with its current simplifying assumptions, which were 
introduced in order to avoid 'over parametrization' and make it predictive and falsifyable, 
will have to be reflned in order to reproduce such data with sufficient accuracy. 

Acknowledgment. We thank an anonymous referee for useful suggestions and T. 
Kruhler for making available to us tabulated data of the optical lightcurves of GRB 071031 
measured with GROND. 
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Table 1. The AG parameters used in the CB model description of the smooth X-ray AG 
of Swift GRBs with superimposed flares. The CB model spectral index obtained from the 
temporal shape of the AG and th at inferred from th e measured XRT spectrum as reported 
in the Swift lightcurve repository ( lEvans et al.ll2009l ) are compared in the last two columns. 
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Tabic 2. The parameters of the ICFs and SRFs used in the CB model description of 
X-ray lightcurves of Swift GRBs. Values of parameters that are not well determined by a 

best fit are inserted within parentheses 
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Ar [sJ 


hj^ykev \l a\s 


GRB050502B 


IGF SRF 


485.7 


308.3 


0.27 keV 


13750 


28950 


6013 s 




SRF 


58950 


18911 


3350 s 








GRB060418 


IGF,IGF 


60 


12 


13keV 


9.9 


118 


0.13keV 


GRB060206 


SRF SRF 


1.19 


32.2 


235 s 


653 


3427 


3645 s 




SRF 


609359 


245909 


1980979 s 








SHB050724 


IGF,IGF 


1.58 


114 


1 08 keV 


469 


1422 


1.32 keV 




IGF 


9757 


49512 


3 53 keV 








GRB050916 


SRF IGF 








17836 


1335 


91 keV 




IGF 


9757 


49512 


3.53 keV 








GRB050916 


SRFJGF 








17836 


1335 


91 keV 


GRB060526 


SRF IGF 








233 


42.5 


15 keV 




IGF,SRF 


280 


25.7 


60.2 keV 


85 


106 


288 s 


GRB060929 


SRF,IGF 








472 


45 


117 keV 


GRB071031 


IGF,IGF 


84.8 


46.2 


13.2 keV 


131 


33.8 


0.74 keV 




IGF,IGF 


181.7 


30.5 


4.3 keV 


228 


40.6 


14 63 keV 




IGF,IGF 


328 


176 


2.74 keV 


572 


173.6 


19 74 keV 




IGF 


4351 


1201 


4.5 keV 










SRF, SRF 


16743 


5463 


4.2 s 


33241 


(53238) 


<C a 


GRB070704 


SRF IGF 


f63) 


(0.35) 


f25 


472 


45.22 


117 keV 


GRB080506 


IGF SRF 


130 


36.7 


14.4 kcV 


211 


43 


171 s 




IGF,IGF 


433 


f41) 


fl59 keV) 


684 s 


fl60) 


fl56 keV) 


GRB080607 


IGS,IGF 


65.9 


8.6 


145 keV 


113.9 


(15.4) 


(350 keV) 




SRF,SRF 


222 


332 


0.6 s 


652 


1018 


335 s 


GRB081203A 


SRF 


79 


10.8 


(0) 








GRB080810 


IGF,SRF, 


74 


16 


25.5 keV 


98.6 


4.6 


34 s 




IGF,SRF 


196 


(14) 


(203 keV) 


265 


208 


(7.2 s) 


GRB080906 


SRF,SRF, 


(49) 


(1) 


(47.7 s) 


157 


27.5 


10.75 s 




SRF,IGF 


490 


8.58 


200 s 


836 


78 


145 keV 


GRB081102 


SRF,IGF 


63 


11 


42.7 s 


914 


88 


111 keV 


GRB090417B 


SRF,IGF 


336.2 


203 


82.6 s 


1227.8 


431 


1.96 keV 


GRB090618 


IGF,SRF 


65.9 


12.9 




(0) 


884 


(0) 



- 28 - 



Table 2 — Continued 
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Table 3. The parameters of the CB model description of the early-time optical 
SRFs/peaks in GRBs 080319B, 060206, 071003, 071031 and 080330 



Flare to [s] a [s] t^xp [s] /3o 
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Table 4. CB model parameters of the ICS 7 peaks in GRB 990123 
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Table 5. The parameters used in the CB model description of unresolved prompt emission 

SRF in the optical hghtcurves of Swift GRBs. 
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Fig. 1. — Schematic illustration, not in scale, of the typical environment encountered by a 
highly relativistic jet ejected in core collapse SN that escapes from the star formation region 
into the galactic halo. 
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Fig. 2. — C omparison betwee n 0.3-10 keV lightcurves of Swift GRBs (Swift/XRT lightcurve 
repository (lEvans et al.ll2009l )) with early-time X-ray flares and their CB model description 
with the parameters lis ted in Table 1 . Top left (a): Comparison between the Swift 'canon- 
ical X-ray lightcurve' (iNousek et al.l 120061 ) of GRB 060729 and its CB model description. 
Top right (b): Zoom on the prompt emission and fast decay phase in a. The CB model 
lightcurve of the prompt emission consists of a sum of 5 ICFs. Middle left (c): Com- 
parison between the evolution of the effective photon spectral index of GRB 060729 in the 
0.3-10 keV X-ray band as inferred from the Swift XRT observations by IZhang et al.l (120071 ) 
and that inferred from the CB model lightcurve. Middle right (d): Comparison between 
the non-canonical X-ray lightcurve of GRB 060418 and its CB model description. Bottom 
left (e): Comparison between the XRT and CB model lightcurves of SHB 050724. Bot- 
tom right (f): Comparison between the evolution of t he eff ective photon spectral index in 
the 0.3-10 keV X-ray band as inferred by IZhang et al.l (120071 ) from the observations of SHB 
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Fig. 3. — Comparison between the 0.3-10 KeV X-ray lightcurves of GRBs with prominent 
flares in their X-ray afterglow that were measured by the Swift XR T and reported in the 



Swift /XRT lightcurve repository http://www.swift.ac.uk/xri_curves/ (lEvans et al.ll2009l ) and 
their CB model description with the parameters listed in Table 1. Top left (a): GRB 
050502B. Top right (b): Comparison between the evolution of the effective photon sp ectral 
index of GRB 050502B in the 0.3-10 keV X-ray band as inferred by (jZhang et al.ll2007l ) from 
the Swift XRT observations and that inferred from the CB model lightcurve. Note the 
similarity between this Fig. and Figs. 2c,f. Middle left (c): GRB 060526. Middle right 

(d) : Zoom on the prominent X-ray flare in the early-time AG of GRB 060526. Bottom left 

(e) : GRB 070704. Bottom right (f): Zoom on the prominent X-ray flare in the early-time 
AG of GRB 070704. 
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Fig. 4. — Comparison between the 0.3-10 KeV X-ray lightcurves of GRBs with prominent 
flares in their X-ray afterglow that were measured by the Swift XR T and reported in the 
Swift /XRT lightcurve repository http://www.swift.ac.uk/xrt_curves/ (lEvans et al.ll2009l ) and 
their CB model description with the parameters listed in Table 1. Top left (a): GRB 
080506. Top right (b): Zoom on the prominent X-ray flare in the early-time AG of GRB 
080506. Middle left (c): GRB 080607. Middle right (d): Zoom on the prominent X-ray 
flare in the early-time AG of GRB 080607. Bottom left (e): GRB 080810. Bottom right 
(f): Zoom on the prominent X-ray flare in the early-time AG of GRB 080810. 
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Fig. 5. — Comparison between the 0.3-10 KeV X-ray lightcurves of GRBs with prominent 
flares in their X-ray afterglow th at were measured by the Swift XRT and reported in the 
Swift/XRT lightcurve repository (lEvans et al.ll2009l ) and their CB model description with 
the parameters listed in Table 1. Top left (a): GRB 080906. Top right (b): Zoom on 
the prominent X-ray flare in the early-time AG of GRB 080906. Middle left (c): GRB 
090417B. Middle right (d): Zoom on the prominent X-ray flare in the early-time AG 
of GRB 090417B. Bottom left (e): GRB 090621A. Bottom right (f): Zoom on the 
prominent X-ray flare in the early-time AG of GRB 09062 lA. 
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Fig. 6. — Comparison between the 0.3-10 KeV X-ray lightcurves of GRBs with prominent 
flares in their X-ray afterglow th at were measured by the Swift XRT and reported in the 
Swift/XRT lightcurve repository (lEvans et al.ll2009l ) and their CB model description with 
the parameters listed in Table 1. Top left (a): GRB 060929. Top right (b): Zoom on the 
flare in the AG of GRB 060929. Middle left (c): GRB 081102. Middle right (d): GRB 
050916. Bottom left (e): GRB 090709A. Bottom right (f): Zoom on the early-time 
X-ray flares in the AG of GRB 090709A. 
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Fig. 7. — Comparison between the observed early-time X-ray and optical lightcurves of very 
bright GRBs and their CB model descriptions. GRB 080319B - Top left (a): 0.3-10 keV 
XRT lightcurve. Top right (b): R-band lightcurve. Middle left (c): The early-time 
optical flares. Middle right (d): The spectral index lightcurve. GRB 071003 - Bottom 
left (e): The 0.3-10 keV XRT lightcurve. Bottom right (f): The entire R-band lightcurve. 
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Fig. 8. — Comparison between the observed early-time X-ray and optical lightcurves of XRFs 
and their CB model description. Unlike the X-ray flares, the early-time optical flares are 
barely resolvable to separate optical flares. XRF 080330 - Top left (a): The 0.3-10 keV 
XRT lightcurve. Top right (b): The optical (white) lightcurve. XRF 071031- Middle left 
(c): The 0.3-10 keV XRT lightcurve. Middle right (d): Zoom on the early time XRT 
lightcurve. Bottom left (e): The e arly-time optical (white) lightcurve. Bottom right 
(f): Evolution of the optical spectral index. 
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Fig. 9. — Comparison between observed and CB model description of light c urves. GRB 



9901 23 - Top left (a): Comparison between the 20-50 keV BATSE hghtcurve flBriggs et al. 
19991 ) and its CB model description, Eq. ([1]), in terms of 9 ICS peaks + a constant background 
of 3850 counts s~^. Top right (b): Comparison between the V band lightcurve of GRB 
990123 and its CB model description assuming a single CB moving in circumstellar density 
profile cx 1/r^ overtaken by a constant ISM density around an observer time t = 1000 s. The 
prompt flare was not resolved into separate flares. GRB 061007 - Middle left (c): The 
0.3-10 keV XRT lightcurve. Middle right (d): The optical lightcurve with evidence for at 
least two early time overlapping flares. Bottom left (e): The optical (white) lightcurve of 
GRB 030418 - Bottom right (f): The R-band hghtcurve of GRB 090102 with evidence for 
a tail of a prompt optical flare. 
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Fig. 10. — Comparison between the observed early-time X-ray and optical lightcurve of 
GRBs and their CB model description. The X-ray and optical emissions do not appear to 
be correlated. GRB 090618 - Top left (a): The 0.3-10 keV lightcurve. Top right (b): 
The optical R-band lightcurve. GRB 090812 - Middle left (c): The 0.3-10 keV XRT 
lightcurve. Middle right (d): The optical R-band lightcurve with evidence for at least 
two early time overlapping flares. GRB 091029 - Bottom left (e): The 0.3-10 keV X-ray 
lightcurve. Bottom right (f): The R-band lightcurve with possible evidence for unresolved 
optical flares. 
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Fig. 11. — Comparison between the observed X-ray and optical lightcurves of XRFs and 
their CB model description. The early-time optical lightcurves are approximated by a single 
(unresolved) SRF. GRB 050820A - Top left (a): The 0.3-10 keV XRT lightcurve. Top 
right (b): The optical R-band lightcurve. GRB 081203A Middle left (c): The 0.3-10 keV 
XRT lightcurve. Middle right (d): The optical U-band lightcurve. GRB 060206 Bottom 
left (e): The 0.3-10 keV lightcurve with poorly resolved late-time flares. Bottom right 
(f): The R-band lightcurve with well resolved late-time SRFs. 



